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Observation of Kinetic Plasma Jets 
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Under certain conditions an intense kinetic plasma jet is observed to emerge from the apex of 
laboratory simulations of coronal plasma loops. Analytic and numerical models show that these jets 
result from a particle orbit instability in a helical magnetic field whereby magnetic forces radially 
eject rather than confine ions with sufficiently large counter-current axial velocity. 
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ny laD and space plasmas (e.g., solar coronal 
loops 3L spheromaks 0], tokamaks Q], and magnetic 
clouds |4[) are presumed to be magnetic flux tubes filled 
with plasma confined via magnetohydrodynamic (MHD) 
forces. However, confinement can be significantly de- 
graded in ways not predicted by MHD; e.g., in tokamaks, 
ions resulting from neutral beams injected against the 
toroidal current direction (counter-injection) exhibit se- 
vere orbit losses compared to co- injection @, @, 0). Re- 
lated confinement degradation may be the cause of small 
solar corona jets (eg,-, surges) associated with canceling 
magnetic features [g, l£j, [l0( and of coronal streamers em- 



anating from magnetic neutral lines . This Letter re- 
ports that in certain circumstances an ion injected along 
the axis of a magnetic flux tube will be magnetically 
ejected from the flux tube instead of being magnetically 
confined, i.e., the ion will have radially unstable motion 
(RUM). This instability explains the severe orbit losses 
of counter-injected ions in tokamaks and is likely relevant 
to similar situations occurring in the solar corona. The 
instability, labeled as 'Kinetic Plasma Jet' in Fig.[IJ was 
discovered experimentally and then modeled. 

We first outline the physical basis for RUM. Con- 
sider a particle injected with velocity v z o near the axis 
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FIG. 1: (Color online) Kinetic plasma jet emanating from 
an argon laboratory loop. Dashed arch corresponds to initial 
plasma loop at t « 1.0 fj,s as in Fig. SJb). Arrows repre- 
sent lines of sight (los#) used for spectroscopy; arrow widths 
represent ~ 6 mm diameter of lines of sight. 
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FIG. 2: (a) Flux tube geometry used in the model, (b) Ex- 
perimental configuration; B outside the flux loop is weakest 
on the right hand side (indicated as weak field region). 



of a cylindrical flux tube having helical magnetic field 
B =£>00 + B z z. The flux tube geometry is sketched in 
Fig. [U^a) and corresponds to a straightened-out model of 
Fig.[2fb), our laboratory configuration simulating a coro- 
nal loop [12]. The r component of mdv/dt = qv x B is 
mr — mrrp 2 = qr<f>B z — qv z B$. If v z B^ = 0, then ra- 
dial force balance r = gives <f> — —qB z /m = —lu c , i.e., 
the conventional cyclotron orbit. However, if v z B^ ^ 0, 
then radial force-balance r = requires cp 2 + (f>qB z /m — 
v z qB c /,/mr = so no real <p solutions exist if 



Bi 



AmVzB^/qr < 0. 



(1) 



Thus, large negative mv z B^,jqr causes the radially out- 
ward force —qVzB^ to overwhelm the radially inward 
force qr<f>B z so centrifugal force mrqT is unbalanced. 

We next use Hamiltonian arguments to show that 
satisfying Eq. (fTJ) leads to the particle being radially 
ejected from the flux tube, i.e., RUM. In order to model 
the simplest nontrivial situation, both B z and the ax- 
ial current density J z are assumed uniform within the 
flux tube so in the flux tube the vector potential is 
A(r) = cj>B z r/2 — z[i J z r 2 /4, the axial flux is <I> = B z -Kr 2 , 



and the axial current is I = J, 



Using the Lagrangian 



L = m(v 2 + r 2 cp 2 + v 2 )/2 + qr<\>A$ + qv z A z , the canonical 
momenta = dL/d<fi and P z = dL/dv z are 



mr 2 4> + qr 2 B z /2, P z = mv z - n qJ z r 2 /4. (2) 



Because <f> and z are ignorable, both P^ and P z are in- 
variants. A particle injected with velocity v z q along the 
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the flux tube axis (i.e., at r = 0) thus has the invariants 

(3) 



^ = 0, P z 



mv z0 . 



Combining Eqs. (0) and (J3]) gives <j) — —ui c /2 and v z — 
v zO +w c ? ,2 A/4 where A = \i§J z jB z = [1qI/<& is related 
to twist. The Hamiltonian H = m(v 2 + r 2 <j) 2 + v 2 ) / 2 can 
be expressed as H = mv 2 /2 + f(r) where 
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(4) 



is an effective potential. On defining £ = |A|r, the di- 
mensionless effective potential U(£) = 2X 2 f(r)/mu> 2 can 
be written as 



ne) = e 2 /4+(A«,oM+c 2 /4) : 



(5) 



Equation (O gives dV/d£ = £/2 + Xv z0 (,/uj c + ^ /A. Near 
the flux tube axis the £ 3 /4 term is negligible, so negative 
<9Vy<9£ near the axis corresponds to having 



S = Xv zQ /lo c + 0.5 < 0. 



(6) 



Our main result is that if S < so dV/d^ < near 
the axis, a particle at r = is on an effective potential 
hill as shown in the S = —0.5 curve in Fig. [3£a) and 
will fall radially out of the flux tube, i.e., RUM. Since 
A = 23^/3^, Eq. is identical to Eq. All mag- 
netic flux tubes have uniform B z and J z near the axis so 
particles with S < will always be on a potential hill and 
experience RUM. Equation ^ can be written in terms 
of experimental parameters as 



S = KIv z0 /q + 0.5 < 



(7) 



where K = iiQirmr 2 j '$ 2 is positive, showing that RUM 
requires Iv z o/q < 0, i.e., counter-current flow. Because 
K (x m, ions have a much lower threshold for RUM than 
electrons. 

Figure G^a) plots U(£) given by Eq. © for S = -0.5 
and 1.5, while Fig.[3jb) plots trajectories calculated from 

direct numerical integration of mr =qvx (^B^cj) + B z 2j 

for a particle with S = —0.5 (i.e., v z0 — —lo c /X) starting 
at the down arrow and for a particle with S = 1.5 (i.e., 
v z q = +uj c /X) starting at the up arrow. Injection at 
Xx = Xy = 10~ 6 is used so a particle does not start 
exactly at the top of a potential hill. To approximate 
the weak field to the right of the flux loop sketched in 
Fig. [2](b) , an exponentially decaying B z in the current- 
free external region is used in the numerical calculation. 
Figure [3](b) shows that the £ = —0.5 particle is ejected 
from the flux tube (i.e., RUM), whereas the S = 1.5 
particle remains on the flux tube axis (i.e., is confined). 

Our experimental configuration [T3 |. sketched in Fig. 
0b), involves top and bottom electrodes (respectively 
cathode and anode) mounted on the end dome of a large 
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FIG. 3: (a) S — 1.5 gives 'valley' (i.e., stable) effective po- 
tential U(£) while S — —0.5 gives 'hill' (i.e., unstable), (b) 
Numerically calculated particle trajectories, (c) Probe mea- 
surement of flux tube magnetic field showing that field is weak 
on right side as sketched in Fig. HJb). 



vacuum chamber (base pressure ~ 10~ 7 mbar). The ex- 
perimental sequence is: (i) slow (~ 10 ms) electromag- 
nets behind the electrodes create an initial arched vac- 
uum magnetic field, (ii) a fast (~ 1 ms) gas valve in- 
jects neutral gas from orifices in the electrodes, (iii) a 
~ 1 kJ, 59 /jF capacitor switched across the electrodes 
breaks down the neutral gas, (iv) a bright plasma loop 
appears. The 10-20 /js dynamical evolution of this loop 
is imaged [l3[ by a fast digital framing camera. Detailed 
measurements in a similar experiment 14| showed that 



the bulk plasma in the flux loop is many orders of mag- 
nitude denser than the injected pre-breakdown neutral 
gas and results from fast MHD ingestion into the loop of 
orifice-originating plasma [TBI . Figure [3fc) shows mag- 
netic probe measurements of flux tube B^ and B z 
as functions of distance L from the flux tube axis in 
the direction away from the electrode plane (data de- 
convolved as in Ref. Q); the magnetic field amplitude 
decays rapidly to the right (corresponding to weak field 
region in Fig. [2fb)). 

Figure [4] shows the evolution of Ar plasma loops for 
different injected gas mass M n and different flux $ as 
a function of time measured from breakdown. M n was 
determined using a thermocouple gauge and has only a 
relative meaning because the plasma shot, being much 
shorter than the gas puff time, uses only a fraction of 
M n . Since the plasma has very low impedance, the ca- 
pacitor acts approximately as a current source. This is 
consistent with the observation that I and hence X/uj c are 
essentially unaffected when M n is varied. However, the 
plasma velocity v z q is observed to be strongly dependent 
on M n with higher v z q observed at smaller M n . 

Figure |4ja) corresponds to M n = 4.9 mg, $ = 0.75 
mWb; Fig. gfb) to M n = 2.3 mg, $ = 0.75 mWb; and 
Fig. g(c) to M n = 2.3 mg, $ = 0.25 mWb. In the 
first two frames of Figs. [4ja-c) the plasma has a smooth 
arch shape; I is low at this stage and the plasma follows 
the half-torus profile of the initial vacuum magnetic field 
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FIG. 4: (Color online) Evolution of laboratory plasma loops 
at (a) high M„, (b) low M„, and (c) low $. Movies placed in 
Ref. |13l ] show the evolutions even more dramatically. 



spanning the electrodes. Then, as / increases, the plasma 
minor radius decreases due to self-pinching while the ma- 
jor radius increases due to the hoop force [17| associated 
with the poloidal magnetic field produced by /. While 
this is happening, the loop undergoes MHD kink insta- 
bility and the projection of the writhed loop axis results 
in a cusp-like dip at the apex 



17j . In the second frame 



(i.e., 2.5 /xs) of Figs. |4{a,b) a finger- like stream of plasma 
emerges near the top (i.e., near cathode) of the loops. In 
Fig. HJb) , which corresponds to low M n and hence high 
v z q , the stream moves toward the ground plane near the 
cathode and leads to a major disruption in /. As also 
seen in Fig. IDJb) , this is followed by the detachment of 
the loop from the electrodes and, for t > 4 /us, formation 
of a plasma jet propagating far to the right of the elec- 
trodes (see also Fig. [1| into the weak field region (i.e., 
to the right in Figs. EJb) and[3^c)). A significant drop 
in / is observed during the detachment phase as well as 
an associated upward voltage spike. From this time on, 
/ commutates to a new shorter path between the elec- 
trodes, while the detached plasma jet propagates away 
from the electrodes. When <I> is lowered as shown in Fig. 
BJc), two critical stages of the detachment are clearly 
seen in the 4.0-7.0 /xs frames, specifically the loop first 
detaches from the cathode and then from the anode to 
form an intense plasma jet. Figures HJa-c) also display 
S estimated using measured cathode region quantities in 
Eq. (O at 2.5 /its (i.e., just before detachment) and indi- 
cate that the plasma jet development in Figs. [4jb, c) is 
associated with having S < 0. 

Equation (|7|) shows that only ions with v z q being large 
and negative relative to / can have S < 0. Measure- 
ments (discussed below) indicate that near the cathode 
ions with large negative v z q (~ 40-60 km/s) indeed exist. 
The slowing-down time (> 100 fis) of these fast ions by 
the plasma (density ~ 10 20 m~ 3 ) is much longer than the 
plasma duration, therefore collisions cannot affect their 
orbits. Since the ion contribution to electric current nec- 
essarily flows in the same direction as the current, ion 
drift motion associated with electric current cannot ac- 



FIG. 5: Measured spectra of an Ar + line (rest-frame wave- 
length Ao = 434.806 nm shown by vertical lines). Velocity is 
v — c(A — Ao)/Ao where A is measurement wavelength, c is 
speed of light. Lines of sight (los #) are shown in Fig. [T] (a) 
Spectra from the cathode (los #1) and anode (los #2) re- 
gions for t — 0-0.5 fis. (b) Spectra from the kinetic jet region 
along los #3, #4, and #5 for t = 2.5-5.5 its, 7.0-9.0 fis, and 
6.0-18.0 fis respectively. 



count for the observed large negative v z q. Furthermore, 
because the measured \v z q\ greatly exceeds the Ar + ther- 
mal speed vt ~2-5 km/s estimated using the spectro- 
scopically determined Tj ~1-10 eV, neither can ion ther- 
mal motion account for the observed large negative v z q. 
However, there does exist a mechanism capable of ac- 
celerating ions to high velocities either parallel or anti- 
parallel to /. This mechanism 14. 
gradients of B^ provide an MHD force ~ 
accelerates plasma from regions of large B\ to regions of 
small £??, i.e., acceleration occurs from both foot-points 



shows that axial 
- -dBl/dz that 



of a flux loop towards the apex if the flux loop minor 
radius is smaller at the foot-points than at the apex (see 



15(|) . The resulting velocity is 



detailed discussion in Ref. 
v z ~ B$l \fJiQffwii, consistent with higher ion axial ve- 
locity observed at smaller neutral gas injection pressures. 



Ar + Doppler velocity measurements have been made 
using aim monochromator with a gated intensified CCD 
camera with fiber/lens coupling system. The spectra dis- 
played in Fig. [5] show velocity components along lines of 
sight (los) indicated in Fig.fJJby "los #". The los #1 and 
#2 spectra in Fig. [5^a) show that both cathode and an- 
ode emission lines are blue-shifted, confirming suprather- 
mal ion flow from both cathode and anode towards the 
apex as predicted by Ref. [HI]. This outflow is seen in 
camera images as 30-60 km/s bright fronts propagating 
away from both electrodes along the flux loop axis to- 
wards its apex [Hj ■ Figure [5ja) shows a large ion veloc- 
ity component (~ 40 km/s ion beam for los #1) mov- 
ing away from the cathode while Fig. [5](b) shows spectra 
measured in the kinetic jet region and indicates a ~ 50 
km/s ion beam for los #3. Because, as seen in Fig. [T] 
los #1 and los #3 make different angles relative to the 
respective flow directions being measured, ion beam ve- 
locities between different los #'s cannot be quantitatively 
compared, i.e., the 50 km/s los #3 ion beam in Fig.[5fb) 
cannot be interpreted as a 10 km/s acceleration of the 40 



4 



0.5 r 
0.0 - 

-1.0 - 

-2.0 - 

-3.0 1 
0.2 





A 


4 






■ 





















-J-it---t- 



0.4 0.6 
$ (mWb) 



0.8 



FIG. 6: Measured 5 for various plasma configurations. 
Filled/open symbols represent high/low M n respectively, 
number of sides in symbols represent capacitor charging volt- 
age in kV, and arrow heads indicate plasma shots where ki- 
netic jets are observed. Existence of kinetic jets has an excel- 
lent correlation with S being negative. 



have much larger orbits (i.e., broader valley-type effec- 
tive potential as in Fig.|3(a)) than co-injected ions. While 
coronal loops are unlikely to have S < due to their small 
A (~ 1CP 8 ni" 1 ) [l a l . je ts associated with canceling mag- 
netic features [Sj, |{J, [lfj and coronal streamers [11| ema- 
nating near magnetic neutral lines are both produced in 
extremely low magnetic field regions where S < could 
occur and RUM may be operative. 

In summary, an instability has been demonstrated 
where ions are magnetically ejected from a flux tube. 
Ejection occurs when ions move opposite to the current 
with a sufficiently large axial velocity. We thank A. H. 
Boozer for pointing out a relationship between RUM and 
neutral beam counter-injection and D. Felt for technical 
assistance. Supported by US DOE and by NSF. 



km/s los #1 ion beam in Fig.[5ja). 

Figure [6] shows the results of a parameter scan of I, 
and v z q performed to determine the S dependence of the 
instability onset. v z q is determined from plasma front 
motion in the camera images [l3J. The S values in Fig. [5] 
were calculated using cathode region parameters in Eq. 
(|7p for a large number of argon plasma loops (r ~ 8 mm) ; 
S = 0.5 is the upper bound for the observed negative 
Iv z o. Kinetic jet instability, shown by arrow-heads in Fig. 
[SI occurs only when S < indicating excellent agreement 
with the RUM onset prediction. 

The plasma loops used for Fig. [6] have already un- 
dergone MHD kink instability [l8j since all have |A| > 
An/L ~ 60 m _1 , where L ~ 0.2 m is the loop length. 
The loops produce kinetic jets only when S < show- 
ing that RUM is a kinetic, rather than MHD, instability. 
The kinetic nature is also evident from the high velocity 
beams in Fig. [5] and from the kinetic jet appearing in the 
weak field (non-MHD) region as sketched in Fig.[2fb). 

The RUM model explains why counter-injected neutral 
beams in tokamaks have severe orbit losses compared to 
co-injected neutral beams @, @, 0- In particular, Fig. 
10 of Ref. [f| showed that an 80 keV counter-injected 
deuterium beam has severe orbit losses in a B — 0.3 T 
tokamak having safety factor q — 1.25 and major radius 
R ~ 1 m. Since A ~ 2/qR ~ 1.6 m _1 , lu cD = 1.4 X 10 7 
s -1 , and the injection velocity is Vi n j = 2.8 x 10 6 m/s, 
it is seen that S coun ter — 0.5 — Xvi n j/u> C D « 0.2 whereas 
S co = 0.5 + Xvinj/uJcD = 0.8; so, counter-injected ions 
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